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Abstract: The numerical hydrodynamic modelling of beat Cepheid behavior has been a longstanding quest in which purely 
radiative models have failed miserably. We find that beat pulsations occur naturally when turbulent convection is accounted 
for in our hydrodynamics codes. The development of a relaxation code and of a Floquet stability analysis greatly facilitates the 
search for and analysis of beat Cepheid models. The conditions for the occurrence of beat behavior can be understood easily 
and at a fundamental level with the help of amplitude equations. Here a discriminant T> arises whose sign decides whether 
single mode or double mode pulsations can occur in a model, and this T> depends only on the values of the nonlinear coupling 
coefficients between the fundamental and the first overtone modes. For radiative models T> is always found to be negative, but 
with sufficiently strong turbulent convection its sign reverses. 
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The Fourier analysis of the observational data of the 
beat Cepheid light curves and radial velocities shows con- 
stant power in two basic frequencies and in their linear 
combinations which indicates that the stars pulsate in 
two modes (or more if resonances are involved). Since 
the beginning of theoretical Cepheid studies in the early 
1960s numerical hydro dynamical attempts at modelling 
the phenomenon of beat pulsation have failed, and beat 
Cepheids have been a bane of stellar pulsation theory. 

In Cepheids energy is carried through the pulsating en- 
velope to the surface by radiation transport as well as by 
turbulent convection (TC). Even though convection can 
transport almost all the energy in the hydrogen partial 
ionization region, this convection is inefficient in the sense 
that it only mildly affects the structure of the envelope. 
It was thus generally thought that convection, while im- 
portant for providing a red edge to the instability strip, 
would play a minor role the appearance of the nonlinear 
pulsation. Purely radiative models did indeed give good 
overall agreement with the observed light and radial ve- 
locities. However, recently it has become increasingly 
clear that there are a number of severe problems with 
radiative models (Buchler 1998), in addition to their in- 
ability to account for beat behavior. 

We have recently implemented in our hydrodynam- 
ics codes a one dimensional model diffusion equation for 
turbulent energy (Yecko, Kollath & Buchler 1998) sim- 
ilar to those advocated by Stellingwerf (1982), Kuhfuss 
(1986), Gehmeyr & Winkler (1992) and Bono & Stelling- 
werf (1994). In contrast to these authors, however, we 
have developed additional tools that allow us to find beat 
behavior without having to rely on very time-consuming 
and sometimes inconclusive hydrodynamic integrations 
to determine if a model undergoes stable, or steady beat 



pulsations. These are (a) a linear stability analysis which 
yields the frequencies and growth rates of all modes, (b) 
a relaxation method (based on the general algorithm of 
Stellingwerf with the modifications of Kovacs & Buchler, 
1987) to obtain nonlinear periodic pulsations (limit cy- 
cles) when they exist, (c) a stability analysis of the limit 
cycles that gives their (Floquet) stability exponents. 

The ID turbulent diffusion equation, and the concomi- 
tant eddy viscosity, the turbulent pressure and the con- 
vective and turbulent fluxes contain (seven) order unity 
parameters that need to be calibrated through a compar- 
ison to observations. In a first paper (Yecko, Kollath & 
Buchler 1998) in which we performed a broad survey of 
the linear properties of TC Cepheid models we found that 
of these the mixing length, the strengths of the convective 
flux and of the eddy viscosity play a dominant role and 
that broad combinations of these three parameters exist 
that give agreement with the observed widths of both the 
fundamental and first overtone instability strips. 

In this Letter we show that the inclusion of TC pro- 
duces pulsating beat Cepheid models that satisfy the ob- 
servational constraints, in particular those of period ra- 
tios, of modal pulsation amplitudes and of their ratios. 
Furthermore the models are very robust with respect to 
the numerical and physical parameters. 

Our discovery of beat Cepheid models has been par- 
tially serendipitous. When we started to investigate the 
nonlinear pulsations of a typical Small Magellanic Cloud 
Cepheid model (M = 4.OM , L^llOOL©, T eff = 5900K, 
Y = 0.73 and Z — 0.004) with the the turbulent convec- 
tive hydrocode, we encountered beat pulsations that ap- 
peared steady. We use the OPAL opacities of Iglesias & 
Rogers (1996) combined with those of Alexander & Fer- 
guson (1994). The values of the TC parameters - for a 
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definition cf. Yecko et al. (1998) - are a c — i 1 o:a = 0.41, 
a p = 0.667, a t — 1, an = 4, a s = 0.75, a v — 1.2. The 
steadiness of these beat pulsations was confirmed when 
several nonlinear hydrodynamics calculations, each initi- 
ated with a different admixture of fundamental and first 
overtone eigenvectors, converged towards the same final 
steady beat pulsational state. This convergence could 
be corroborated when we extracted the slowly varying 
amplitudes with the help of a time-dependent Fourier 
decomposition, and plotted the resulting phase portraits 
(A (t) vs. Ai(t)) that are shown in Fig. 1 where all 
initializations are seen to converge toward a fixed point 
located at A = 0.0104 and A 1 = 0.0200 (These radial 
displacement amplitudes assume the eigenvectors to be 
normalized to unity at the stellar surface, <5r/r„ = l). 








Fig.l: Evolution of the modal amplitudes for different ini- 
tial conditions (marked with a cross). Equal time intervals 
between dots. The open squares denote the unstable fun- 
damental (F), the unstable first overtone (O) and the stable 
double-mode (DM) fixed points. 

While the observed transient behavior of the models 
provides a conclusive proof of the presence of steady beat 
pulsations, it is important to explain and describe the 
behavior on a more fundamental level. The phase por- 
trait of Fig. 1 is very similar to those found for nonres- 
onant mode interaction on the basis of amplitude equa- 
tions. (Buchler & Kovacs 1986, 1987, hereafter BK86 
and BK87) . We show here that indeed the nonlinear be- 
havior of the hydrodynamical model pulsations can be 
understood very simply that way. 

The amplitudes of the two nonresonantly interacting 
modes obey remarkably simple equations 

^ = A (k - qooAl - qoiAj) (la) 

^- = A 1 (ki - q w A 2 - qil A\) (16) 

These amplitude equations are 'normal forms' and are 
therefore generic for any dynamical system in which two 



modes interact nonresonantly. The assumptions underly- 
ing these amplitude equations are satisfied for Cephcids: 
(a) The lowest modes (fundamental and first overtone 
here) are weakly nonadiabatic, i.e. the ratios of linear 
growth rates k to periods are small, a condition that is 
readily confirmed by our linear stability analysis; (b) The 
pulsations are weakly nonlinear which allows a trunca- 
tion of the amplitude equations in the lowest permissible 
(third) order; weak nonlinearity can be established by 
comparing the linear and nonlinear periods which differ 
less than a tenth of a percent. Furthermore both the non- 
linear self-saturation coefficients qoo and qn as well as the 
cross-coupling coefficients qoi and gio have always been 
found to be positive in Ccpheid models so that amplitude 
saturation can occur in third order, and it is sufficient to 
keep terms up to cubic in the amplitudes, (c) In the range 
of interest there is no important low order resonance be- 
tween the fundamental and the first overtone modes, and 
possibly a higher mode. 

In the following discussion we look at the regime where 
both modes are linearly unstable, k > and K\ > 0. 
Eqs. (1) then have two single mode fixed points. The 
amplitude of the single mode fundamental (0) fixed point 
is Ao = \Jk-o/ qoo and its coefficient for stability to first 
overtone perturbations is Ruo) = Ki — <Zio^o • I n our 
notation a positive coefficient implies growth and thus 
instability. The corresponding first overtone (1) limit 
cycle amplitude is A\ = \j K\jq\\ and its coefficient of 
stability to fundamental perturbations is Ko(i) — K o ~ 
qoiA 2 . The k's, when multiplied by the periods Pk of 
their limit cycles, are equal to the corresponding Floquet 
exponents (Buchler, Moskalik & Kovacs 1991). 

Eqs. (1) can also have a double mode Exed point whose 
amplitudes satisfy A\ DM = R 0(1)qil /V < A%, A\ DM = 
«i(0)9oo/f < A\, where V = q oqn - qoiqio- This fixed 
point exists provided A^ DM > and A\ DM > 0. Then, 
if V < 0, the double mode limit cycle is unstable. Sta- 
ble pulsations occur either in the fundamental or first 
overtone, and the pulsational mode is determined by the 
evolutionary history of the model (hysteresis). On the 
other hand, if V > 0, the double mode fixed point is sta- 
ble, and steady double mode pulsations occur (BK86). 
One can show that these conditions are equivalent to re- 
quiring Ko(i) > and Ki(o) > 0, conditions which also im- 
ply that both single mode limit cycles (fundamental and 
first overtone) are individually unstable. This validates 
Stcllingwerf 's (1975) suggestion that the simultaneous in- 
stability of the fundamental and the first overtone leads 
to steady beat pulsations (in the absence of a resonance). 
It provides an economical tool to search for double mode 
behavior, because we can now relatively easily compute 
single mode limit cycles and their stability. 

As a further confirmation that the nonresonant sce- 
nario applies to the pulsating Cepheid model, we have 
determined the coefficients of Eqs. (1) as in BK87 by fit- 
ting time-dependent solutions of these equations to the 
temporal variation of the amplitudes in their approach to 
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the limit cycle as shown in Fig. 1. The fitted trajectories 
in the phase portrait are practically undistinguishablc 
from the hydro results, confirming the applicability and 
accuracy of the amplitude equation formalism and the 
absence of any relevant resonances. 

We mention in passing that the expression 'dou- 
ble mode Cepheids' is often used cavalierly for beat 
Cepheids. Since no additional, resonant overtone is in- 
volved in the beat pulsations, the latter are thus truly 
double mode pulsations. 

With the relaxation code we are able to compute both 
the fundamental and the first overtone limit cycles with 
their respective amplitudes and Floquet stability expo- 
nents Ai(o) = fo^i(o) an d A (!) = Pi/t (i)- The above dis- 
cussion then shows that from these four quantities we can 
extract the four nonlinear qjk coefficients when we have 
already computed the linear periods and growth rates. 
The values we obtain this way for this beat Cepheid 
model agree quite well with those that we obtain from 
the fit described in the previous paragraph. Note that 
these two determinations rely on independent numerical 
hydrodynamical input, the first on two periodic limit cy- 
cles (that are linearly unstable), the second on transient 
evolution toward the stable double mode pulsation. 

In order to investigate the robustness of the observed 
beat behavior we now explore the pulsational behavior 
of a sequence of Cepheid models in which the effective 
temperature of the equilibrium modes of the sequence 
varies from T e jj = 6200K to 5800 K. Such a sequence is 
approximately along an evolutionary path. The eddy vis- 
cosity parameter a v is treated as an additional variable 
parameter to explore the effect of TC on the behavior. 
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Fig. 2: The limit cycle stability coefficients {dT l ) along a 
Cepheid sequence (i.e. as a function of T e //) for various eddy 
viscosity strengths a v ; fundamental limit cycles /tuo) (full 
circles), first overtone limit cycles /t (i) (open circles). 

In Fig. 2 the stability coefficients of the sequence are 
plotted versus T e ff, with open/filled circles for those 
of the fundamental/overtone single mode cycles. The 
curves are labelled with the corresponding strengths a v of 
the eddy viscosity. As discussed above we expect double 
mode behavior where both Floquet exponents are posi- 



tive. (The stability exponents due to perturbations with 
other modes are always smaller in this sequence and are 
therefore irrelevant here). For the low value of a„ = 0.5 
(dotted lines) the two stability coefficients are never pos- 
itive simultaneously, thus excluding double-mode behav- 
ior. On the other hand, for a„ = 1.2 a double mode 
region appears between T e //~ 5875 - 5915 K and for 
a„ = 2.0 this broadens to T e// ~ 5965 - 6050 K. 

How does turbulent convection bring about double 
mode behavior? Fig. 2 shows that, in the region of in- 
terest, an increase in the turbulent eddy viscosity causes 
a rapid decrease in the stability of the fundamental limit 
cycle (ki(o), filled circles), but an increase in that of the 
first overtone limit cycle (ko(i)> open circles). This de- 
scription, though, does not tell us whether it is the effect 
of TC on the linear k's or on the nonlinear q's, or on 
both, that is responsible for the beat pulsations. 

Table 1 shows the variation with a v of the relevant 
model quantities, viz. the linear growth rates, the nonlin- 
ear coupling coefficients, the discriminant T> = <7oo<Zn — 
9oi9iOj the amplitude of the fundamental cycle and its 
stability coefficient with respect to overtone perturba- 
tions, and the same for the first overtone. 

The necessary condition for stable double mode pulsa- 
tions, viz. V > 0, is never found to be satisfied in radia- 
tive models. In these models the cross-coupling always 
dominates over the self-saturation coefficients. Table 1 
shows that an increase in the strength of the eddy vis- 
cosity causes goo and qn to increase faster than q i and 
<7io. It is therefore the resultant change in the sign of V 
that makes double mode behavior possible for sufficiently 
large a v . 

The condition for beat behavior is thus seen to be 
rather subtle in that it involves the effects of convection 
beyond the linear regime for which it seems difficult to 
give a 'simple' physical explanation. 

Fig. 3 gives the overall modal selection picture in the 
a v - T e ff plane. The linear edges of the instability region 
(«o = and K\ = 0) are shown as dashed lines. By com- 
puting the fundamental and first overtone limit cycles for 
a number of a„ and T e ff values, by interpolation, we can 
obtain «o(i) or ^i(o) as a function of a v and T e ff, and in 
particular the loci where they vanish. The solid curves 
give the nonlinear pulsation edges and are marked ORE 
and FBE. 

It is straightforward to show that if the two linear 
growth rates vanish at the same point, the four curves 
will intersect in a single point on this diagram, that we 
label critical point. 

The curve marked OBE is the linear blue edge of the 
first overtone mode and it coincides with the overtone 
nonlinear blue edge up to and on the left of the crit- 
ical point. The linear fundamental blue edge becomes 
also the fundamental blue edge above the critical point. 
Above the line ORE we have k (i) >0 and the first over- 
tone limit cycle is unstable. Below the line FBE the 
quantity «i(o) > and the fundamental limit cycle is 
unstable. Thus in the region marked 'dm' both single 
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mode limit cycles are unstable, and this is the region of 
double-mode pulsation. In the small triangular region at 
the bottom, on the other hand, both limit cycles are sta- 
ble, and either fundamental or first overtone limit cycles 
can occur. 
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Fig. 3: Modal selection in the a v - T e ff plane. 

In summary, stable first overtone pulsations occur in 
the dotted region, delineated by the lines OBE and 
ORE. The fundamental limit cycle is stable in the re- 
gion marked by open squares, delineated by FBE and 
FRE (not shown on the far right). This figure makes it 
particularly evident how TC favors double mode pulsa- 
tions and why all efforts with radiative codes have failed 
in modelling beat Cepheids. 

We have seen that when TC effects are sufficiently 
large then the Cepheids should run into the double mode 
regime in both their crossings of the instability strip. Fur- 
thermore, as a Cepheid crosses the double mode regime 
redward, say, the first overtone amplitude should grad- 
ually go to zero while the fundamental amplitude in- 
creases from zero to the value it attains as a fundamental 
mode Cepheid (BK86) . The question arises whether this 
nonresonant scenario that is derived from the amplitude 
equations is in agreement with the observations. 

The four SMC beat Cepheids from the EROS survey 
(analyzed by Beaulicu and reproduced in Buchlcr 1998) 
all have the same amplitude ratios, Aq/Ai ~ 0.45, a pri- 
ori in disagreement with the nonresonant scenario shown 
in Fig. 1. of BK86 that suggests that Cepheids with all 
amplitudes ratios should occur. 

In Fig. 4 we display the behavior of the component 
modal amplitudes of the beat Cepheid models for the 
a v = 1.2 sequence of Fig. 2. The amplitudes of the stable 
single mode limit cycles are shown as solid lines with solid 
dots for the fundamental and open dots for the first over- 
tone, and as dashed lines where they are unstable. The 
fundamental and first overtone component amplitudes 
of the stable double mode pulsators are shown as solid 
squares and open diamonds, respectively. It is seen that 



although the modal amplitudes do indeed vary continu- 
ously throughout as the double mode regime is traversed, 
the behavior is very rapid near the cooler side. The rea- 
son for this unexpected behavior is that the q's are not 
constant in this sequence, and what is more, they vary 
in such a way that V happens go through zero around 
5850 K. It is the presence of this nearby pole that causes 
a change in the curvature of A . 
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Fig. 4: Sequence of models (same as in top of Fig. 2). Am- 
plitudes of the fundamental and of the first overtone limit 
cycles (lines), solid/dashed where stable/unstable. Double 
mode component amplitudes of the double mode (fundamen- 
tal: solid squares and first overtone: open diamonds). 



According to Fig. 4 it is therefore much more likely to 
find beat Cepheids in the slowly varying regime where the 
ratio Aq/Ai & 0.5. The computed behavior of the modal 
amplitudes is thus in agreement with the observed SMC 
Cepheids, and the nonresonant scenario is consistent with 
the observations. 

We have demonstrated that turbulent convection leads 
naturally to beat behavior in Cepheids, which does not 
occur with purely radiative models. The reason is that 
the nonlinear effects of TC dissipation can create a re- 
gion in which both the fundamental and the first over- 
tone cycles are unstable, and the model undergoes stable 
double mode pulsations. At a more basic level the ampli- 
tude equation formalism shows that turbulent convection 
modifies the nonlinear coupling between the fundamental 
and first overtone modes in such as way as to allow beat 
behavior. 

The development of a relaxation code (TC) to find 
periodic pulsations, and a Floquet stability analysis of 
these limit cycles has made the search quite efficient, and 
a broader survey of beat Cepheids, with wide ranges of 
metallicities is in progress. This will also search for beat 
Cepheid models that pulsate in the first and second over- 
tones. 
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(T-026031), AKP (96/2-412 2,1) and by the NSF (AST95- 
28338) at UF. Two of us (JPB) and (ZK) thank the French 
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TABLE I.- Nonlinear Coupling Coefficients (5900 K model). 
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